Application of two dimensional Prenkel-Kontorova model to nanotribology 
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A two-dimensional Frenkel-Kontorova model is set up. Its application to the tribology is consid- 
ered. The materials and the commensurability between two layers strongly affect the static friction 
force. It is found that the static friction force is larger between two layer of same materials than 
that for different materials. For two-dimensional case the averaged static friction force is larger for 
the uncommensurate case than that for the commensurate case, which is completely different from 
one-dimensional case. The directions of the propagation of the center of mass and the external 
driving force are usually different except at some special symmetric directions. The possibility to 
obtain superlubricity is suggested. 
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Driven dynamics of a system of interacting atoms is an 
interesting physical problem. It has important applica- 
tions in mass and charge transport phenomena in solids 
and crystal surfaces. One of important applications has 
emerged in tribology studies, where a thin atomic layer 
is confined between two substrates which move respect 
to one anotherUQ. Fundamental understanding of fric- 
tion is vitally important in many areas of science and 
technology, ranging from nanotribology to crack prop- 
agation of earthquake dynamics [3]. Progress has also 
been made on how to tune the intrinsic frictional forces 
between a sliding subject and the underlying substrate, 
through interface modification down to the molecular or 
atom scale[3, 0]- Such important understanding at the 
microscopic level, in turn, is expected to serve as impor- 
tant guidance in the design of smart materials with de- 
sirable lubricant properties for industrial and biomedical 
applications [l, Q . 

The application of driven Frenkel-Kontorova (FK) type 
model has received an increasing interest as a possible 
interpretative tool to understand the complex field of 
nanotribology [1]. One dimensional FK model has been 
extensively studied recently [1, However, in friction 
systems, an importance of higher dimensionality has been 
emphasized, which makes the superlubricity (the state of 
vanishing friction) appear much more easily TH 12, T^ . 
It is found in the experimental results that by measuring 
the atomic-scale friction force as a function of rotational 
angle between two contacting layers superlubricity may 
appear in certain misfit angle between two lavers[lllll2|. 
Friction forces completely vanish when, for example, in- 
commensurately contacting surfaces sliding against each 
other. In such contact, the ratio between the lattice units 
of the surfaces is irrational along the sliding direction, so 
each individual atom receives different amounts of force 
from different directions. These forces consequently off- 
set each other, resulting in zero friction. This offsetting 
of forces is made possible by the continuous motion of 
atoms, which is the basic principle behind superlubricity. 



Higher dimensionality of system is, therefore, crucial for 
atoms to move continuously. It is, therefore, necessary 
to extend the ID FK model to the higher dimensions to 
obtain a more real model which can be realized in exper- 
iment. 

For this reason, we consider an upper layer in which 
there are N x M atoms and they are arranged on a 
2D square lattice. We first investigate the position 
and the velocity of an arbitrary (n, m)th atom, where 
n = 1, 2, • • • , A^, TO = 1, 2, • • • , M. Its position can be 
expressed by r„^m = (a^n.m, where a;„,™ and ?/„,™ 

are the positions in x and y directions respectively. For 
this arbitrary (n, TO)th atom, the interactions among the 
nearest and the next nearest neighbors are considered 
among the upper layers. The interatomic interaction 
potential is chosen to be of the simple harmonic form 

Vint = Y.r,j -f [(a;ij-a;„,m-/)^-l-(2/ij -y„,m-0^], with a 
strength K , a natural equilibrium spacing I = a between 
the nearest neighbor atoms of (ri -I- 1, TO)th, (n — 1, TO)th, 
(n, m -|- l)th, (n, to — l)th, and I — \f2a between the next 
nearest neighbor atoms of TO-|-l)th, (n— 1, TO-|-l)th, 

(n -|- 1, TO — l)th, [n— 1, TO — l)th for the upper layer. 

However, this arbitrary (n, m)th atom not only inter- 
acts with each other among the particles in the upper 
layer, but also with the lower layer through a 2D periodic 
substrate potential which depends on the lattice struc- 
ture in the lower layer. For different materials we choose 
different substrate potential. For generality, we choose 
two kinds of substrate potentials. One is with a square 
lattice symmetry, Vsqare = ^[K) - cos - cos 
where / is the magnitude of the adhesive force be- 
tween the two layers, Vb is a constant and the length 
h is the natural equilibrium spacing of the lower layer. 
The other is with a hexagonal symmetry, Vhexagonai = 
I(y„_cos22^cos^). 

For general case of the system, in which the orienta- 
tions of the two layers do not match, we rotate the two 
layers with respect to each other by an arbitrary misfit 
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angle 9. Then 



The po- 



a; \ _ /cosf —sinU\ x 
y) ^ \sin0 cos9 J \yy 
sition of this arbitrary (n, m)th atom Tn m satisfies the 
foUowing equation of motion, 



int 



(1) 



where 7 is a phenomenology viscous damping coeffi- 
cient. It can be thought of representing degrees of free- 
dom in real physical systems which are not explicitly in- 
cluded in our model (e.g., vibrational or electronic ex- 
citations in the substrate). We use dimensionless vari- 
ables, where identical mass m = 1 is assigned to each 
atom. The frequencies of atomic vibrations are isotropic, 
i.e., LUx = (jjy = •\/27r/, the characteristic time scale is 
To = ^27r//. Fext = {Fext COS a, FextsiT^ a) is the exter- 
nal driving force, and a is the angle between directions 
of Fext and the unit vector of x axis. 

A fourth-order Runge-Kutta algorithm was imple- 
mented to solve Eq.(l) numerically. The averaged ve- 
locity is defined as v = nxai ^\'=!i 'j=i ^ < where 
<> denotes the time average. We modeled N x M atoms 
of upper layer and choose the periodic boundary condi- 
tion to enforce a fixed density condition for the system. 
XM+i — xi + Ma, UN+i = yi + Na. The initial condi- 
tion we choose here is that the velocity of each atom is 
zero and the position of each atom is at its equilibrium 
position. For our system we let iV = 12, M = 12. In 
order to study how the static friction force varies with 
the different materials of the lower layer, the numerical 
simulation of the mobility as a function of the driving 
force and the static friction force as a function of differ- 
ent system parameters of both upper and lower layers are 
presented. 
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FIG. 1: The mobility < v > as a function of driving force 
Fe^t for / = 1, a = 6 = 1, m = 1, is: = 1, 7 = 0.7, 6 = 40° and 
different values of a = 0°,25°,45°. 
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FIG. 2: The numerical results of Fs as a function of a in 
different 0. (a)The commensurate case with the square sym- 
metry substrate potential. (b)The golden mean case with the 
square symmetry substrate potential, (c) The spiral mean 
case with the square symmetry substrate potential, (d) The 
commensurate case with the hexagonal symmetry substrate 
potential, (e) The golden mean case with the hexagonal sym- 
metry substrate potential, (f) The spiral mean case with the 
hexagonal symmetry substrate potential. 



We define a new parameter /3 = b/a^ H EB]. For 
simplicity, we study three typical cases of /3 = 1 (com- 



mensurate case), P = (golden mean case), and 

/3 = 0.755 (spiral mean case), respectively, for both 
square and hexagonal symmetric substrate potential. 

Fig.l show the numerical results of the average chain 
velocity as a function of the driving force F^xt for the 
case of = 40°, X = 1, / = 1, and 7 = 0.7 with 
different values of a. It is for the commensurate case 
with square symmetric substrate potential. We note from 
Fig.l that the average velocity v is zero if the external 
force is less than Fg (the static friction force). As the 
force Fext increases adiabatically, the system undergoes 
a sharp transition from the pinned phase to the running 
crystal phase. Our numerical results indicate that for 
other values of 9 the system also transfers pinning state 
directly to the sliding state, although Fg is different for 
different 9. Similar results are obtained for the case of 
hexagonal symmetric substrate potential. It seems that 
the magnitude of Fs depends on both parameters of a 
and 9. Therefore, the static friction force depends on the 
external driving force direction and the misfit angle. We 
also find from our numerical results that Fs decreases as 
either / decreases or K increases. 

In order to understand the anisotropic characters of 
the system the numerical results of Fs as a function of 
a in different 9 are given in Fig. 2. It is noted that, for 
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the case of square symmetric substrate potential with 
P = 1, Fs strongly depends on the parameter of 6, espe- 
cially at 6* = the static friction force Fg is much larger 
than that for other values of 9. Fg also depend on the 
parameter of a. But for the golden mean case and the 
spiral mean case the static friction forces Fg depend on 
but not as much as that of the commensurate case of 
/3 = 1 (see Figs.2(a)-2(c)). The variations of averaged 
friction forces with respect to 9 are not as much as that 
of commensurate case. 

Fig. 2 also show how the static friction force varies with 
respect to different materials of the substrate. We choose 
two kinds of substrate materials. One is with square 
symmetry and the other is with hexagonal symmetry. 
We find that Fg varies with different materials of the 
lower layer. In order to know how the static friction 
force depends on the different materials of the lower layer 
we give the averaged values of Fg for different substrate 
potentials and different parameters of (3 which is shown 
by Table 1. 



Table 1 Averaged static friction forces 
obtained from the numerical values of Fig. 2. 



Upper layer 


Lower layer 


/3 


Averaged Fg 






1 


0.521 




Square lattice 


0.755 


0.766 


Square lattice 




0.618 


1.058 






1 


0.46 




Hexagonal lattice 


0.755 


0.713 






0.618 


1.013 



It seems from Table 1 that the static friction force Fg 
is larger between two layers of same materials than that 
for different materials which is in agreements with exper- 
imental results [TtI. 

Meanwhile, we note from Fig. 2 that the averaged fric- 
tion forces for two typical uncommensurate cases are 
larger than that of the commensurate case of /3 = 1 for 
both square and hexagonal symmetric substrate poten- 
tials. It means that for 2D case the averaged friction 
force is larger for the uncommensurate case than that for 
the commensurate case. This result is completely differ- 
ent from ID case, in which case the friction force is larger 
for the commensurate case than that for the uncommen- 
surate case [13, [l^ . 

Now we analyze the atom trajectories for the case of 
/ 1, 7 ^ 0.7, a = 1, 6 = 1, /s: 1, 61 = 40°, and a = 25° 
with the square symmetric substrate potential shown in 
Fig.l. As the external driving force increases, each atom 
of the system may move from its equilibrium position. 
Fig. 3(a) present the trajectories of the (5,4)th, (5,5)th 
and (5,6)th atoms at Fext = 0.1. It seems that there 
are small displacements around their equilibrium posi- 
tions for these atoms. However, the displacements are 
much less than lattice spacing a. Each atom move near 



its equilibrium position. However, the averaged atom ve- 
locity is still approximately zero. This state corresponds 
to a pinned state. When the driving force increases, the 
state will be transferred into a sliding state of moving 
crystal, see Fig. 3(b) for Fgxt — 1-82, and Fig. 3(c) for 
Fext — 1-84, respectively. Fig. 3(b) show that the direc- 
tion of average atom velocity is not same as that of the 
external driving force. The propagation direction of each 
atom is not a constant, but there are vibrations around 
its average direction. It actually presents a atom mo- 
tion in a solitonlike fashion. We observe a formation of 
moving kinks for each atom. It suggest that for differ- 
ent driving forces the directions of average atom velocity 
are different. Fig. 3(c) show a complete crystalline state 
in which case the directions between driving force and 
the average atom velocity are same. Although there are 
small vibrations between atom trajectory and the direc- 
tion of the driving force, it becomes smaller and smaller 
as the driving force increases. 
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FIG. 3: The trajectories of selected atoms for different exter- 
nal driving forces, (a) The trajectories of the atoms of (5, 4)th, 
(5, 5)th and (5, 6)th at Fext ~ 0.1. (b) The trajectories of the 
atoms of (5, l)th, (5, 2)th, (5, 3)th, (5,4)th, (5, 5)th, (5,6)th, 
(5, 7)th, (5, 8)th, (5, 9)th, and (5, 10)th at F^^t = 1.82, (c) The 
trajectories of the atoms of (5, l)th, (5, 2)th, (5, 3)th, (5, 4)th, 
(5,5)th, (5,6)th, (5,7)th, (5,8)th, (5,9)th, and (5, 10)th at 
Fext = 1.84. 

Fig. 3 only show the trajectories of particular atoms of 
the system. In order to know the trajectory of mass cen- 
ter (cm.) of the whole system, the differences between 
the direction of external driving force and the direction 
of trajectory of the center of mass are given in Fig. 4. We 
define a parameter of a' to represent the angle between 
directions of trajectory of cm. and the unit vector of x 
axis. The dependence of a' on the external driving force 
are shown for 9 = 0°,45°,40° in Fig.4(a), 4(b), 4(c), re- 
spectively. It is found that the propagation direction of 
the center of mass is same as that of the external driving 
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force when a = 45° for both Q = and B = 45° due to 
the symmetry of the system. However, two directions are 
different for other cases, even for the cases of = 0° and 
Q = 45°. For other values of d two directions are differ- 
ent for any values of a. We also note that if the external 
driving force is below the static friction force, although 
it is not in sliding case, the propagation direction of the 
center of mass is always in the direction of o! = 45°. 



Wigner crystal(20j. metallic dots 21, 23|, magnetic bub- 
ble arrays[23j, etc. 




Driving force (F ) 



FIG. 4: The angle ol between the direction of trajectory of 
the center of mass and unit vector of x axis as a function of 
the magnitude of the external driving force F^^^t. (&)Q = 0°. 
(b)e = 45°. (c)6l = 40°. 

In a conclusion, we find that the averaged static fric- 
tion force is larger between two layers of same materials 
than that for different materials. For 2D case the av- 
eraged friction force is larger for uncommensurate case 
than that for commensurate case. To obtain superlu- 
bricity, we may choose different materials, but with com- 
mensurate ratio between two layers with larger stiffness 
strength K of upper layer and smaller magnitude of ad- 
hesive force between two layers /, which may be a case 
of real system such as diamond-like carbon filmQ, about 
which superlubricity may be realized. The directions of 
the propagation of the center of mass and the external 
driving force are usually different. We may devise a ex- 
periment to verify this result in the future and it may 
be useful to the many fields of condensed physics, such 
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